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ABSTRACT 

Upcoming ground- and space-based observatories such as the DES, the LSST, the JDEM 
concepts and the SKA, promise to dramatically increase the size of strong gravitational lens 
samples. A significant fraction of the systems are expected to be time delay lenses. Due to the 
sensitivity of time delays to cosmological parameters, this increase in sample size opens up 
new avenues of parameter constraint. Many of the existing lensing degeneracies, e.g. between 
the Hubble parameter Hq and the slope of the density profile of the lens, become less of an 
issue with large samples since the distributions of a number of parameters are predictable, and 
can be incorporated into an analysis, thus helping to lessen the degeneracy. Assuming a mean 
galaxy density profile that does not evolve with redshift, a ACDM cosmology, and Gaussian 
distributions for bulk parameters describing the lens and source populations, we generate 
synthetic lens catalogues and examine the relationship between constraints on the - ^^A 
plane and Hq with increasing lens sample size. We find that, with sample sizes of ~ 400 time 
delay lenses, useful constraints can be obtained for and JIa with approximately similar 
levels of precision as from the best of other methods. In addition, sample sizes of '^lOO time 
delay systems yield estimates of Hq with errors of only a couple of percent, exceeding the 
level of precision from current best estimates such as the HST Key Project. As such, strong 
lensing promises to become one of the premier methods for constraining the Hubble constant 
in the future. We note that insufficient prior knowledge of the lens samples employed in the 
analysis, via under or overestimates in the mean values of the sample distributions, results in 
broadening of constraints. This highlights the need for sound prior knowledge of the sample 
before useful cosmological constraints can be obtained from large time delay samples. 
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1 INTRODUCTION 

Strong gravitational lensing provides us with a unique and com- 
pelling probe of cosmology. In particular, it has been known for 
some years that time delays between multiple images of certain 
strong len sing systems a llows for constraints on Ho, the Hubble 
constant, ( lRefsdallll964l) . due to the dependence of the time de- 
lay scaling on this parameter. The constraints on Hq by means 
of strong lensing time delays has historically had mixed suc- 
cess, mainly due to its degen eracy with the matter density pro- 
file of the lensing galaxy (e.g . Falco, Gorenstein & Shapir 

1 



iGorenstein, Shapiro & Fa"icoiri98^ ; ICorless. Dobke & Kindf2008h 
However, despite a number of remaining systematic uncertainties. 
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we h ave recently witnessed increasi ng precision from this method 
(e.g. ISaha et aflHoO^ : lOguril l2007h together with general agree- 
ment with the HST K ey Project value of 72 ± 8 kms"^ Mpc~^ 
jpreedman et alj200lh . 

Lens systems and their time delays, however, are also sensitive 
to other cosmological parameters, specifically the matter density of 
the Universe, Qm, the dark energy density Qoe, the dark energy 
equation of state parameter, w, and the curvature Qk- Our current 
best u nderstanding of the universe is that w = -1 (e.g. Spergel et al] 
l2007h . consistent with dark energy as a cosmological constant; 
from hereon we refer to JIde as f^A. Together with Hq, these pa- 
rameters enter lens models through the angular diameter distances 
between the source, lens, and observer, as discussed below. 

Previous analyses of lens statistics have already yielded useful 
constraints on Ha- A combination of the Cosmic Lens All-Sky Sur- 
vey (CLASS; lirowne et alj[2003l ; iMvers et al.ll2003l) with recent 
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Sloan Digital Sky Survey (SDSS) data on the local velocity disper- 
sion distribution function of E/SO galaxies teemardi et"5]|2003h 
was used to provide a constraint of Q,a ~ 0.74 — 0.78 (assum- 
ing flatness), and a 95% confidence upper bound of Q.a < 0.86, in 
good agreement with results from the redshift-mag nitude relation 
of Ty pe la supemovae and from CMB experiments jMitchell et al.l 
l2005h . 

The sensitivity of the time delays is at its strongest when 
considering the Hubble constant, allowing useful constrai nts to be 
obtained even when analysing only one system (s ee e.g. ICourbinI 
l2003l : iKochanek & Schechteiil2004i ; |jacksodl2007l for a review of 



Ho estimates from lensing). This is not the case for the remaining 
cosmological parameters however, since the sensitivity of the time 
delays to Qm, ^a, and is relatively weak. This results in the 
normal practice of fixing these parameters to the concordance cos- 
mology as derived from a combination of measurements (e.g. SNla 
[Tonry et al. 200^; CMB|Spergel et al. 2007), and focusing on using 
lenses to constrain Ho- 

The situation could quickly change in the near future. The past 
two decades have witnessed a golden age of astronomy and cosmol- 
ogy, with observatories such as the Hubble Space Telescope (HST), 
Chandra, Keck, the Very Large Telescope (VLT), Very Large Array 
(VLA), and the Wilkinson Microwave Anisotropy Probe (WMAP), 
to name but a few, providing us with discoveries that have changed 
our perception of the cosmos from the smallest to the largest scales. 
However, we are now beginning to enter a new era, as these great 
observatories and instruments are either replaced or augmented 
by a next generation of space- and ground-based telescopes. In- 
cluded amongst these are the Large Synoptic Survey Telescope 
(LSST;|lyezic 2007), GAIA erryman 2001), Dark UNiverse Ex- 
plorer (DUNE; Refregier et al. 2006), the European Extremely 
Large Telescope (E-ELT ; Gilmozzi & Spyromilio,2007 ), the Thirty 
Meter Telescope (TMT; ISilva et al.ll2007l). the Jo int Dark Energy 
Mission (e.g. DESTINY teenford & Lauej |2006h . ADEPT Q and 
SNAP jjel insky 2007) concepts), and the Square Kilometre Array 
(SKA: ICarilli & Rawlingsl 120041) . This range of instruments, and 
others not mentioned above, will have a direct impact on the cos- 
mology that can be done with strong lensing. Predictions for the 
number of galaxy-scale strong lensing candidates vary between the 
instruments, but conservative estimates put the combined tally into 
the many tens of thousands. One particular example, the proposed 
NASA/DOE Joint Dark Energy Mission concept SNAP, is pre- 
dicted to discover bet ween 5,000-50,000 lensed galaxies , and 100- 
1,000 lensed quasars ( iMarshall. Blandford & S^l2005h . a signif- 
icant fraction of the latter having corresponding time delays (cur- 
rently ~ 15% of known lens systems are time delay lenses). Some 
telescopes will make discoveries of lens systems fortuitously, in 
much the same way as is being done with analysis of the Sloan 
Digital Sky Survey (SDSS) data, while others will have dedicated 
programs searching for lens candidates - notably the LSST with 
its high cadence observations providing time delay monitoring on a 
sub-weekly basis. Instruments such as the Square Kilometre Array, 
which lie a littl e further in the future, could also yie ld thousands of 
lens systems ( Koopmans . Browne & Jacksodl2004h . The proposed 
30 — lOOm-class ELTs and the TMT will enable the lens galaxies 
themselves to be studied in exquisite detail, helping to constrain 
their velocity dispersions and density distributions. We highlight 
that samples provided by the above instruments, even when consid- 
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ered purely on an individual basis, will dwarf the current samples 
of ~ 20 time delay lenses. 

This prompts the question as to how such an increased sam- 
ple of time delays will be used to benefit cosmological studies, in 
particular in the constraint of cosmological parameters such as fim, 
SI A and Ho- Indeed, while it is true that the angular diameter dis- 
tances are only weakly sensitive to the cosmological parameters 
Sim, S^A, and Qk, a vastly increased sample size would compen- 
sate for this factor. Additionally, many of the existing lensing de- 
generacies, e.g. between Ho and the slope of the density profile 
of the lens, should become less of an issue with large samples 
of lenses since many parameters should tend towards a distribu- 
tion that can be incorporated into an analysis. This would provide 
another independent constraint on those cosmological parameters, 
and in combination with other methods, using CMB measurements, 
SNla, and X-ray clusters, would result in further increases in pre- 
cision. Hence, the aim of this work is to investigate the relationship 
between constraints on the S7m - S7a plane and on Ho as a function 
of increasing lens time delay sample size. 

This paper is organised as follows. We begin in §2 by describ- 
ing the generation of large synthetic data sets and justifying the ex- 
act distributions of parameters used. §3 will describe how we then 
take these data sets and analyse them to obtain meaningful con- 
straints on the cosmological parameters discussed. §4 presents the 
key findings from our investigation, highlighting constraints in the 
S^m ~ S^A plane and briefly examining the impact that increased 
sample size can also have on the determination of Ho- Finally, §5 
will provide a discussion of our results and present conclusions. 



2 A SYNTHETIC LENS SAMPLE 

We begin with the generation of a synthetic data set upon which 
an analysis can be performed. To accomplish this it is important to 
put into place a justifiable prescription for predicted features of fu- 
ture lens samples, specifically what characteristic value parameters 
that define the sample may take. For simplicity, below and in the 
analysis that follows, we assume Gaussian distributions for quanti- 
ties describing the lens and source population characteristics (lens 
and source redshifts, lens galaxy ellipticity, velocity dispersion and 
slope of surface mass density profile, and convergence and shear of 
lens galaxy environment). In addition, we do not account for co- 
variance between parameters. Clearly, while Gaussian distributions 
and neglect of covariance serves as a fairly good approximation to 
some of these quantities, for others this is a very crude represen- 
tation. The implications for these approximations are discussed in 
§5. 

The key parameters that will govern our sample are as follows: 
1) Lens and source redshift (zi, Zs); We currently have good con- 
straints on source qua sar redshifts from surveys such as SDSS (e.g. 
ISchneideretar]|2003b and the Two-Degree Field QSG Redshift 



IScmieider et al.l l/UU^I) ana tne Iwo-Ueg ree tieia IjSU Keasmtt 
Survey (2QZ) (e.g. Tang & Zhand 12003) . which include sample 
sizes in excess of 40,000 QSOs. Similarly, since we know that the 
vast majority of lens galaxies are early-types, we ca n again make 
use o f large studies such as the DEEP survey (e.g IWeiner et al.l 
I2OO5I) and the ESO-Sculptor survey (e.g. lde Lapparent et alj2003h . 
both of which consist of ~ 500 early-type galaxies (note that 
SDSS is too shallow for this purpose). These samples indicate a 
peak in the lens and source redshifts of around zi = 0.4-0.6 and 
Zs = 1.70. However, we also need to take into consideration the 
greater sensitivity of future observatories. For example, predictions 
for the peak lens and source redshifts from both the SNAP mis- 
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Figure 1. Esti mates of pred icted parameter distributions for future lens samples. Gaussian properties from t op left to bottom right : (a) Lens redsliift, = 0.5, 
a = 0.15, (e.g |Weineretalj [2005; de Lapparent et al. 2003 D ; (b) Source redshi ft, = 2.00, a = 0.75, (e.g. lSchneider et al.ll2003t iTang & Zhandl2005l) : (c) 
ElUpticity, = 0.25, c = 0.1, <Bak & Statleri200Q : iKhairul A lam & Rvden'2002!) ; (d) Velocity dispersion, = 230 km s"!, cr = 50 km s ^ (e.g. lSheldon et all 
12004) : (e) Density profile slope, = 2.00, a = 0.50, (e.g. [Rusin, Kochanek & Keeton 2003; Koopmans et al. 2 0oi) ; (f) Group shear, = 0.075, a = 0-025 , 
(e.g. fpalal & Watso j2004lHolder & Schechtej|2003l) ; (g) Group convergence, ^l = 0.075, cr = 0.025, (e.g. .Dalar& Watsorj2004lHolder & Schechtej|2003l) 



sion concept and LSST [f| are clo ser to zi ~ 0.5 and Zs ~ 2.0 
( [Marshall. Blandford & SakdbOOSi) . In addition, the existing sam- 
ple of lenses have lens and source redshift distributions more in line 
with these slightly higher estimates 0. The input redshift distribu- 
tions for both lenses and sources are estimated from the above data, 
with widths of ct = 0. 15 and a = 0.75 respectively, and are shown 
in Fig.[TJi andflj). 

2) Galaxy ellipticity (e); As stated, the majority of lens galaxies ex- 
hibit some form of ellipticity. We input an estimate o f the predicted 
distrib utio ns in ellipticity on the sky as t aken from iBak & Statlen 
( l2000h and lKhairul Alam & Rvdei] J2002h . where they measure the 
distributions of almost 1800 elliptical galaxies giving a mean of e 
^ 0.25 (as defined via the axis ratio of the light; e = 1 - b/a). This 
is a very similar level of ellipticity as found in the Sloan Lens ACS 
(SLACS) survey sample (e ~ 0.24) as measured in both the mass 
via strong lens model ling and in the l ight, although this is from a 
sample size of just 15 jTreu et al.l2006h . The input distribution used 
in the generated data set, with corresponding width of ct = 0.1, is 
displayed in Fig.[T};. 

3) Velocity dispersion ( ct„ j; Once again, we use data from the SDSS 
survey taking velocity dispersion distributions from a sample of 
almost 50,000 early-type elliptical galaxies ( Sheldon et al. 20(3)- 
Here they split the available sample in half, giving velocity disper- 
sion distributions for both low and high ct„ ellipticals, the means of 
which were ct„ = 143 kms~^ and ct„ = 215 kms~^ respectively. 
Lenses will have a bias towards the high mass systems since the 
lensing cross-section scales with greater mass/velocity dispersion. 
Hence the higher ct„ sample represents a better estimate for the 
distribution of lens galaxy velocity dispersions. Indeed, this is also 
in keeping with the SLACS sample of 15 lenses which had a peak 



^ See|http://www.ls st.org | 

^ See castles lens database; |http;//www.cfa.harvard.edu/castles/] 



velocity dispersion of between ~ 240-280 km s~ . The greater sen- 
sitivity of future instruments may however enable the detection of 
lensing by slightly lower mass systems. Hence we take an approx- 
imate mid-point between the SDSS and SLACS sample and use a 
mean of 230 km s^^ in our synthetic sample. The form of this dis- 
tribution, with corresponding width of ct = 50 km s~^, is shown in 

4) Density profile slope (rj); While often a contentious is- 
sue, it seems that most lens galaxies studied so far do in- 
deed have a nearly isothermal ma ss distribution. Studies suc h as 
iRusin. Kochanek & Keetoi] ( |2003|) and lKoopmans et ahl ( l2006h use 
lensing data (combined with stellar dynamics data in the latter 
work) from 22 and 15 lens systems respectively to obtain val- 
ues very near ?7 = 2 for a density distribution that follows a 
form ptot oc poT^^. Indeed, it is likely that on average the 
slope will converge in a large sample (i.e ^ 100) to a value 
which, on current indications, is around rj = 2. There are how- 
ever lenses that do appear to have density profil es that are signif- 
icantly steeper or shallower than isotherinal (e. g. ^Kochanekj |2002 ; 
iTreu & KooDmansll2002l ; lKochanek et alJl2006l ;IAuger et all boOj ; 
Read, Saha & Ma ccio 2007) for a number of possible reasons (e.g. 
iDobke King""&Fe llhauer 2007; Auger 2008; Read et al. 2007). In 
order to accommodate such lenses we introduce a significant spread 
to the distribution of ct = 0.50 around our isothermal mean of r; = 
2. This is shown in Fig.[T^. We will discuss the factors effecting the 
density profile and its degeneracies further in §2.1. 

5) Group convergence and shear (k, 7); It is harder to estimate 
the average contributions from the group environment, particularly 
from observational data. For this study we have chosen to con- 
sider expected distrib utions as calculated from N-body simulations 
of group galaxies. In lOalal & Watsoril ( l2004h , both the group con- 
vergence and shear are fou nd at levels of ~ 0.03 - 0.05, whereas 
iHolder & Schechted ( I2OO3I) find external shear values closer to ~ 
0.11. Again, lenses may be biased toward systems with a higher 
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group convergence and hence may be more likely to be found in 
future surveys (e.g . iMomcheva et alj[2006l : iFassnacht et alJllOOd : 
lAuger et"aDl2007an . Similarly, shear values may be correspond- 
ingly higher - some modelling examples of the time delay lens 
PGl 115+080 require a minimum of 7 = 0. 05 to successfully re- 
produce image positions l lWitt&Maolll997h . Given these factors, 
for this study we have chosen to increase the mean group conver- 
gence over those found in the N-body investigations to k = 0.075, 
and to select the mid-point of the shear from the two studies, i.e. 
7 = 0.075. We display our input distributions in Figs.[Tf and[T^ 
respectively. Both have corresponding widths of a = 0.025. 

We note that all the distributions given above and shown in 
Fig.[T]are Gaussian fits using the available surveys and studies, and 
allow us to generate sample data. While it is not expected that sim- 
ilar Gaussian fits to future samples will change significantly from 
the distributions given, variation is inevitable. To help compensate 
for this effect, we have generally increased the spread of the distri- 
butions over that taken from the literature, as highlighted earlier for 
the case of the density profile slope. 

Once we have obtained reasonable estimates of parameter dis- 
tributions, we can begin constructing synthetic data sets. Unlike the 
distributions discussed above, the source position /32) for each 
system is selected from a uniform distribution. For the purpose of 
this investigation, the data sets will consist of the following ob- 
servables - image positions and the corresponding time delays, as 
generated from a particular set of input parameters; velocity dis- 
persion ((T„), density profile slope (rf), ellipticity (e), source posi- 
tion j32), source and lens redshift {zs, zi), the Hubble constant 
(Ho), external shear and convergence (7, k). Additional constraints 
can be obtained with the inclusion of fluxes but traditionally it has 
been hard to make models fit well to the observations (e.g. due to 
micro-lensing, substructure effects etc.), so we exclude flux ratio 
constraints. We use an elliptical power-law potential to generate 
time delay surfaces, then locate image positions and calculate the 
corresponding time delays. The potential has the form 



[{1 ~ e)el + {1 + e)e. 
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where 
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with Ds and Dis being the angular diameter distances between 
the observer and source, and the lens and source respectively. 
The terms 61 and O2 are the angular image positions in the im- 
age plane. We use an elliptical power-law potential model for 
speed and simplicity and for the reason that when investigat- 
ing low ellipticity galaxies as is done in this study, the surface 
mass density corresponding to a power law potential model is 
virtually indistinguishable from that of a true surface mass den- 
sity model (jKassiola & Kovner 1993). We model the contribu- 
tion of the external group convergence and shear as follows (e.g. 



iKeeton & Zabludofd[2004l) 



[k + ■ycos2{e - e-.f)] , 



(3) 



where the positions 9i and 62 are related to the position angle, 6, 
via 9i = rcosd and 82 = rsinO with r = 9f + 9^, and is the 
shear angle. 



2.1 Galaxy density profiles and associated degeneracies 

As highlighted above, it is difficult to simultaneously constrain 
the density profile slope (and its evolution) and Hq, due to their 
strong degeneracy. We mentioned in our discussions above that cer- 
tain time delay lenses appear to have steeper density profiles than 
those determined from the SLAGS sample. In principle this could 
result from a systematic slope offset between objects with well- 
measured slopes where time delay measurements are not possible, 
and lenses with time delays but poorly determined slopes. We note, 
however, there are also time delay systems which yield Ho > 72 
km s^^ Mpc^^ when modelled with an isothe rmal lens, imply - 
ing a shallower density profile e.g. B0218+357 (York et al. 2005k 
81608-^656 iKo opmans et al. 2003), HE l 104-1805 (Ofek&Maoj 
|2003|) . Q0142-100 jAkhunov et aLll2008l) . In addition, there is not 
any direct evidence that the SLAGS sample is insensitive to such 
systematic errors, and therefore the SLAGS distribution of may 
represent the true underlying distribution including any systemati c 
effects. When considering the analysis performed in lOguril ( l2007h . 
we also note that on an individual basis, before any combination 
of data was performed, the Ho values are scattered in an approxi- 
mately symmetric way around Ho = 72 km s^^ Mpc^^, and that 
this could be attributed to the approximately symmetric nature of 
the systematic errors that effect mass distribution profiles. Indeed, 
based on the best currently available empirical data, it seems that a 
broad Gaussian distribution based around a value of 77 = 2 is a fair 
starting point for the analysis that follows. Another factor is that 
the measured profile slope could possibly vary with radius. This 
would be a concern if the slope is calculated by comparing lens- 
ing c onstraints to dynamical information (e.g., lTreu & Koopmand 
I2OO4I) . The velocity dispersion measurements are made at R < 
Reff, while the lensing information comes from the Einstein ring 
radius, Reiu- An effective slope, rjefj, calculated between these 
radii may not match the local value of t) in the annulus between the 
lensed images if there is a significant difference between Rgf f and 
Rsin - Derivations of the Hubble Constant are particularly sensitive 
to any variation in the local value of rj, and to achieve percent level 
accuracy one would need any slope correlation to be approximately 
dr;/dR < 0.02 kpc^^ around the region of the Einstein radius. How- 
ever, it appears from the SLAGS study at least, that there are no 
correlations between the density slope and the effective radius, the 
Einstein radius, the distance between th e two radii, and the mean 
of the two radii dKoopmans et al.ir2006h . These measurements are 
based on 15 lenses in the SLAGS sample, and an analysis based 
on a larger sample of ~100 systems will help to rule out any such 
correlation. 

We note that there is also a possibility that galaxy density pro- 
files could evolve with redshift/time, and could therefore be cor- 
related with certain factors e.g. with the angular diameter distance 
to the lens. The best observational evidence on evolution comes 
from co mbining the SLAGS and LSD (Lenses Structure and Dy- 
namics; iKoopmans & Treull2003l) data sets; assuming a simple lin- 
ear dependence of 77 on z gives d77/dz = 0.23±0. 16 (Ict) over the 
range z = 0.08 — 1.01, marginally consistent with no evolution. 
However, the mean of th e Lenses Structure and Dynamics (LSD; 
IKoopmans & Trej|2003l) sample is outside the joint-analysis range 
from the SLAGS lenses, which could either be a redshift effect, 
since the LSD lenses are at a higher redshift range than those of 
the SLAGS study, or simply be due to the small number statistics 
involved in the LSD study. To beat down the ±0. 16 error to a level 
of around ±0.02 would require approximately 80 additional lens 
systems. More important however is the redshift coverage of that 
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sample. Due to the limited depth o f SPS S much of the SLACS 
sample used in the iKoopmans et alj ( l2006h study is only out to a 
redshift of z 0.3. A deeper sample of lenses, as will be discov- 
ered in the next generation of lens surveys, together with an even 
coverage over the z = 0-1 range, will help us understand the true 
evolution of slope vs. redshift. 

The SLACS study now has around 80 lens systems and dy- 
namical mass estimates of higher redshift systems are being de- 
termined - this wi ll help in better understanding the evolu tion of 
the density profile teolton et al.l2008al : lBoiton et al-llOOSbh . Addi- 
tionally, analysis of large N-body simulations may help to establish 
whether any subtle evolution should be expected. Regardless, the 
evolution certainly does not appear to be extreme and for this study 
we assume that there is no correlation between redshift and den- 
sity slope. In the future, with large increases in lens survey data, 
stronger constraints on evolution will be derived and incorporated 
into analyses. 

The steepness of the lens density profile is not the only fac- 
tor that can be degenerate with the Hubble constant. Shape de- 
generacies can also exist, for examp l e in th e form of varying and 
twisting ellipticities (see ISaha et si\ ^200 6*) for discussions). Re- 
cent studies have shown that the projection of intrinsically triaxial 
lens galaxies carmot effect the derivation of the Hubble constant 
jCorless, Dobke & King|[20o3) . Further studies are required in or- 
der to explore other possible shape degeneracies. For now, however, 
it is not entirely clear that the SLACS lenses would not have fairly 
sampled this distribution and at present the best empirical distribu- 
tion of lens galaxy density slopes is approximately Gaussian around 
rj = 2. 

Given larger upcoming samples, it may well be feasible to de- 
fine a pre-selection process for systems that are to be included when 
deriving constraints on cosmological parameters. For example, sys- 
tems in environments with a richness above a certain threshold, 
knowingly p rone to profiles that deviate from isothermality (e.g. 
lAugeill2008b could be excluded. It is possible that these selection 
criteria could introduce other biases into the sample. On the other 
hand, should no measurable differences in results be found between 
the inclusion and exclusion of such systems, their relative impor- 
tance could be better understood. 



3 MODEL AND ANALYTIC METHOD 
3.1 The angular diameter distances 

The dependance on cosmology in this study comes via the 
observer-lens, observer-source and lens-source angular diameter 
distance expressions, Di, Ds and Dia respectively. Angular diam- 
eter distances are related to line-of-sight comoving distances 



dz' 

Ho J, W)' 



(4) 



where the function E(z') depends on Qm, ^a, and ^Ik in the fol- 
lowing way 



E{z) = ^/Qm{l + Z)-' + niiil + ZY + S^A . 



(5) 



Di, Ds, Dis depend on Dc in different ways. When Hk = 0, 
as will be assumed in th e following analysis, consistent with obser- 
vations of the CMB e.g. lSpergel et alJoOOTh . we have 



Diiz) 



Dcjzi) 
l + zi 



(6) 



Dsiz) = 



l + Zs 



Du{z) = — \D,{zs) - Dc{zi) 

1+ Zs 



(7) 
(8) 



3.2 Elliptical power-law lens model 

We model the data set described in §2 with an elliptical power-law 
(EPL) model, the form of the potential being that of expression ([TJ. 
The corresponding two dimensional lens equation is given by the 
lens equation j3 = 6 ~ \7tjj{9), where ^p{9) is the combination of 
the lens potential and the potential contribution from the external 
convergence and shear, i.e. tptot = tpepi + '>Pk,t Hence for the two 
components of /3 we have 



/3i = (l-«: + 7)6li 
'b{l~e)e,{[{l 



and 



[{i-e)ei + {i + e)eiYn 



/32 = (1 - «: - 7)^2 

" 6(1 + - e)el + (1 + e)el]^'^ /bf-^ 



, (9) 



(10) 



The time delay at a particular image position (with respect 
to the light travel time of the unperturbed ray in the absence of 
lensing) calculated by the solution of the lens equations has the 
form (see e.g. Naravan & Bartelmannlll997h 



m 



l+Zl 



DiDs 



o-p\ -i>tot{e) 



(11) 



The resultant time delay At between two image positions is 
therefore the corresponding difference in the light-travel time given 
by O- 

The lens equations {9} and dlOt are solved simultaneously us- 
ing the Newton-Raphson method for a n on-linear system o f equa- 
tions. We employ the MNEWT routine jPress et al.|[l992h which 
requires the two lens equation components and the corresponding 
Jacobian matrix. This routine returns the image positions, 6i, that 
satisfy the lens equations for a given set of free parameters; av, rj, e, 
/3i, P2, Ho, 7, K. For a given set of image solutions and parameters 
the model time delays are then calculated using Jilt . The image po- 
sitions (6i) and corresponding time delays (tj) are then compared 
to those from the synthetic data set of §2 using statistics. 



E 



{6i,Mod — Oi,Obs,ey 



+ 



.7 = 1 



(12) 



^ At{Obs) 



where f70(^obs) and o-At(Obs) ars ths errors in the synthetic observed 
positions and time delays which are taken as conservative esti- 
mates of 1% and 10% respectively. We minimise expression l ll2t by 
employing the POWELL minimisation routine (Press et al. 1993) 
which returns the best fit parameters for any particular system in 
an array comprising [av, rj, e, /3i, 132, Ho, 7, k\, for a given choice 
of f2ni and Q,/^, which is explored on a grid. As noted above, we 
assume a flat Universe with Q.k = 0. The corresponding grids 
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(a) 100 systems (b) 200 systems (c) 400 systems 
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Figure 2. (a) The f2m - f^A plane as constrained from a synthetic data set consisting of 100 time delay lens systems. Plotted are the 1, 2, and 3(t contours. 
We see a strong degeneracy along the Qm + Qa = 1 line. Constraints are approximately similar in scale to those from e.g. early CMB experiments, SNla 
measurements. X-ray clusters, but in comparison to more modem day instruments it provides us with no new limits ; (b) Same as left plot but for 200 systems. 
The degeneracy is still present but the contours have shrunk considerably, particularly along the Cm + S^A = 1 direction ; (c) Constraints are now based upon 
data from 400 systems. There is significant closure in contour size when compared with 100 and 200 systems. This level of parameter constraint is similar in 
precision to modern CMB instruments, and it is now possible to recover the input cosmology (Hm = 0.3, Ha = 0.7) to relatively high precision. 



are combined from each system, allowing for the investigation into 
how the constraints scale with the number of systems combined for 
any given survey data set. 



Table 1. Model prior values as used in expression )13t 



Parameter (p) Mean (pp) Variance (cp) 



3.3 Model priors 

As highlighted in §2, there is a great deal of information that we 
already know a priori concerning the properties of early-type lens 
galaxies. In any future analysis of a lens sample we would not nec- 
essarily know the precise values for parameters such as the slope of 
the density profile, or ellipticity, but it is the intention of this work 
to highlight that the bulk properties of the sample are quite well 
defined, and as the sample size increases, the necessity to know the 
precise values for parameters potentially diminishes accordingly. 
More important is the information about the sample as a whole, and 
what prior knowledge we can put into the prescription to help im- 
prove our constraints on cosmological parameters from lens sam- 
ples. To this end we place some additional Gaussian priors into the 
prescription given above l ll2t . These inc lude priors on e.g. the 
Hubble constant from the HST Key Project jpreedman et alj200lh 
when constraining the f2ni ~ Q.A plane. Information on parameters 
such as velocity dispersion are essentially similar to those given in 
§2, however since the exact distribution of the lens sample would 
not be known in reality, but rather only an estimate of the bulk 
properties, we expand those in §2 and use far broader distributions 
as the actual priors. The values used for the mean and variance of 
each prior on Ho, a^, e, rj, k and 7 is displayed in Table 1, and the 
prior takes the form; 



2 _ (P - Mp) 



(13) 



We note that when fitting two image systems, no prior on 7 is 
employed since such systems typically have low levels of external 
shear. 



Ho 


72.0 


16.0 


av 


230 


200 


e 


0.25 


0.25 


V 


2.00 


2.00 


K 


0.075 


0.15 


7 


0.075 


0.15 



4 ANALYSIS AND RESULTS 

Below we employ our analysis method detailed in the previous sec- 
tion with the synthetically generated lens sample. The input cos- 
mology that is used to generate the lens sample is a ACDM model 
with ilm = 0.3, D.A = 0.7, and h = 0.72. 



4.1 Constraints on the fim - f^A plane 

We incorporate the priors detailed in Table 1 and expression l ll3t 
into our function l ll2t . We begin by analysing the Qm - 
plane and examine the relationship between lens sample size and 
parameter constraints with samples of 100, 200, and 400 lens sys- 
tems. The results are shown in Fig.|2] 

To begin with, Fig.|2^ shows constraints based upon the anal- 
ysis of 100 systems. The predominant degeneracy lies along the 
Jim + ^A = 1 direction. Constraints at this level are approximately 
similar in s cale to those from early CMB experim ents (e.g. Max- 
ima, COBE: |jaffa et al.l200ll:IStomDor et alJ2003l). SNlamea sure- 
ments and X-rav clusters jAllen et al.ll2004lTonrv et alj2003l) , but 
in comparison to modem day CMB instruments it provides us with 
no new limits. In Fig.[2j5 we increase the sample size to 200 sys- 
tems and see a notable reduction in the 2, and 3a contours. This 
level of constraint in the fim - Qa plane is similar to that from 
studies involving more mod ern CMB measurements (e.g. WMAP 
- DRh ISpergel et alj|2003h . The first data release of the WMAP 
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instrument yielded the value fim = 0.29 ± 0.07 for WMAP data 
alone, compared to an error of ± 0.1 on the input value of Qm = 
0.30 for Fig.[2j3 above. Similar error comparisons are found for 
though not explicitly stated in DRl. 

To examine the possible limit of this approach we double the 
number of systems again; Fig. [2}; shows constraints based upon a 
sample of 400 lens systems with their corresponding time delays. 
We now see a significant closure of all contours and recover the in- 
put cosmology (ilm = 0.3, SIa = 0.7) with high precision. For this 
sample size the 1, 2, and 3a errors are smaller than t hose produced 
by m odem CMB experiments (e.g. WMAP - DR3; ISDergel et all 
l2007h . When assuming w = -1 and the HST value for Ho, the third 
data release from the WMAP instrument yields values of fim = 
0.241 ± 0.034 and = 0.716 ± 0.055. In comparison. Fig. [2j; 
above yields similar errors of ± 0.05 on the input values of Om 
= 0.30 and Q,a = 0.70. Naturally, future CMB measurement preci- 
sion is also expected to improve, e.g. Planck ( Xia et al. 2007). At 
this point, we note that even if ultimately the precision from this 
method is found to be less than that from CMB experiments, which 
is probable, tight constraints on parameters should be possible, and 
it is the complimentary nature of such a measurement that is valu- 
able (seelLinder (2004) for discussions). 

So far we have assumed that the bulk properties of the data 
sample being analysed are well known i.e. that the mean of the pri- 
ors match those of the synthetic data set distributions of §2. It is 
probable that this will not always be the case and that some prior 
knowledge is lacking with certain distributions. To investigate the 
effects of this we re-run the analysis discussed above, but include 
a ±1(7 shift in the mean of the priors (using the Gaussian widths 
of §2, Figs.[TJ to simulate a systematic under or over estimate in 
one of the bulk sample properties. For the case when one out of 
the six priors in §3.3 had their mean shifted, the recovered con- 
straints remained the same or very similar as for those displayed in 
Fig.|2] When we introduced systematic shifts in the means of two 
out of the six priors, the overall constraint was reduced to that typ- 
ically obtained from half the original number of lens systems, e.g. 
a ±1(T shift in the means of both velocity dispersion and ellipticity 
resulted in constraints obtained for 400 lens systems that were sim- 
ilar in scale to those from 200 systems obtained with no shift in the 
means. If a severe lack of prior knowledge was assumed in the anal- 
ysis, and three or more of the six priors mismatched the synthetic 
sample, or if the shift was ±2(t or more, the constraints no longer 
minimised on the synthetic sample's starting cosmology of Qm = 
0.3, = 0.7. We highlight, that a ±2(7 shift in e.g. the velocity 
dispersion distribution, would imply a severe lack of prior knowl- 
edge of a lens sample and is hard to foresee. When small shifts of 
< ±1(7 where introduced, very little difference was observed in the 
recovered constraints, even when introduced into multiple priors. 

The key investigation in this particular study, i.e. the constraint 
on cosmological parameters as a function of sample size, depended 
weakly on changes in the exact form of the initial starting distri- 
butions used to create the initial synthetic data set. We tested this 
by creating different synthetic data sets after varying the form of 
the distributions discussed in §2 (i.e. varying the quoted jj, and a of 
Fig.[TJ, and repeating the analysis - in all cases the basic relation- 
ship between the constraint on cosmological parameters as a func- 
tion of sample size remained constant. Far more important in the 
ability to obtain constraints on cosmological parameters via time 
delays is the prior knowledge of the bulk properties of that sam- 
ple (whatever form they happen to initially take), in particular the 
prior knowledge on the mean values of the various distributions, as 
discussed above. 



4.2 Implications for Hq 

We have seen that with lens sample sizes of ~ 400 time delay 
systems it is possible to obtain comparatively high precision con- 
straints in the fim ~ plane. As an independent method for the 
study of cosmological parameters this holds the potential to play 
a useful role in the future. Larger sample sizes also have a direct 
impact on the study of the Hubble constant. The current best value 
for the Hubble c onstant is still the HST Key Project value of 72 ± 8 
km s^^ Mpc' jFreedman et al.l200lh . A number of Ho values are 
found from modelling individual lensing time delay sys tems and 
are in the range of about 50-80 km s~^ Mpc~^ ( see e.g. ICourbinI 
l2003l : lKochanek & Sch echtei"2004 ljacksoij|2007h . although when 
considered as an ensem ble goo d agreeme nt with the HST value is 
found (e.g. IOgurill2007i : ISaha etai1l2006l) . The CMB does not di- 
rectly measure Ho, but it is possible to infer a measu rement (based 
on a number of analysis specific assumptions; see ISpergel et al.l 
l2007h of i7o ~ 73 ± 3 km s'^ Mpc"\ Future mcreases m su- 
pernova samples sizes (e.g. from SNAP and DESTINY) may help 
to increase the precision on Ho . In this section, however, we focus 
on the possible improvements on precision that could one day come 
about due to large increases in the number of strong lensing time 
delays. 

As highlighted in § 1 , it is common practice to fix the cosmol- 
ogy to some pre-defined concordance cos mology as derived from a 
combination of me asurements (e.g. SNlal Tonrv et alj|2003l ; CMB 
ISpergel et al]|2007l) . and focus on obtaining estimates for Ho, and 
some success has been had with this approach. One of the main ob- 
stacles to this method is the degeneracy that exists between Ho and 
the slope of the density profile of the lensing galaxy. As discussed 
however, at some particular sample size the average density pro- 
file should tend towards a predictable distribution that can help in 
breaking the degeneracy and from observations this v alue appears 
to be close to isothermal (e.g. iKoopmans et al.ll2006l) . We explore 
the implications of this concept by investigating the variation of 
Ho parameter constraint with increasing sample size, as done pre- 
viously in the Jim - ^a plane. We remove the HST prior, but keep 
the slope of the density profile as 77 = 2 ± 1 and leave the other 
remaining priors as they are. We display results in Fig.[3] 

We begin with a sample size of 25 systems (Fig.[3ji). Notice 
that the constraints placed on Ho in this case are similar to those 
foun d with similar sized samples from existing observational data 
(e.g. lOguril2007l : ISaha et al.l2006l) , the current-day sample consist- 
ing of ~ 20 time delay lenses, which put values at 70 km s~^ 
Mpc~^ (2(7) and 72 tli km s"^ Mpc~^ (2a), respectively. 

We double the sample size to 50 time delay lenses and show 
the results in Fig.[3j3. The contours shrink notably and we are now 
in the realm of precision that exceeds the HST Key Project value 
(72 ± 8 km s~^ Mpc~^). Increasing the sample size to 100 systems 
moves beyond this level of precision, towards ~l-2% uncertainty 
in the recovered value (Fig.[3]:). 

As was done in the case for constraints on Qm - ^a, shifts 
in the priors were introduced when considering constraints on Ho ■ 
We noted similar broadening of the constraints over those shown 
in Fig.[3]when introducing mean shifts in multiple priors, although 
in this case it had less of a negative effect, i.e. the width did not 
increase to the level found with half the number of systems as was 
the case for shifts discussed above for the flm ~ S7a plane. This is 
possibly due to Ho being more sensitive to time delays than S7m and 
f^A, and hence having more of a capacity of overcoming poor prior 
knowledge. Nevertheless, there was still an observable broadening 
of the contours and this should be noted. 
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(a) 25 systems (b) 50 systems (c) 1 00 systems 




60 64 68 72 76 80 84 88 60 64 68 72 76 80 84 88 60 64 68 72 76 80 84 88 

Hq (km s"'' Mpc"'') Hq (km s"'' Mpc'"') Hp (km s'"' Mpc"'') 

Figure 3. (a) The constraints on Hq based upon 25 time delay lens systems. Horizontal lines mark the formal 1, 2, and 3(T limits (dot-dash, dash and dotted 
respe ctively). Sample size and the constraints obtained are similar in magnitude to those obtained from existing time delay samples (~20 systems, c.f. lOguril 
120071) ; (b) The sample size is now doubled to 50 time delay lenses showing the corresponding reduction in width ; (c) Constraints are now based upon 100 
time delay systems, yielding \a en'ors in the region of '^1-2 % uncertainty in Hq. 



5 DISCUSSION AND CONCLUSIONS 

Future observatories are assured of finding a large number of strong 
gravitational lens systems, which have the potential of providing a 
unique and compelling probe of cosmology. There are a number of 
studies that can be performed via this unique phenomena, one of 
which is the well documented use of constraining the value of the 
Hubble constant. Less well explored, however, are the constraints 
that can be put on fim and f^A via lensing time delays, mainly 
due to the low numbers of such systems available for study at the 
present time. Lens samples promise to grow rapidly once upcoming 
surveys get underway including the DES, LSST, JDEM, DUNE, 
and the SKA, which will uncover many thousands of strong lens 
systems. Investigations into what cosmological information these 
lens samples can yield could help shape the particular survey strate- 
gies themselves. 

Lens statistics have already pr ovided us with usef ul con- 
straints in the Om - plane (e.g. [Mitchell et al.ll200^ . While 
some ~ 100 multiply imaged quasars and radio sources have 
been discovered, a statistical analysis requires a sample from a 
survey that is complete and has well-characterised, homogeneous 
selection criteria. T h e largest such sam ple comes from CLASS 
( lBrowneetani2003l : iMvers et alj|2003h . with a total of 22 con- 
firmed lenses. Future surveys, such as those performed by the LSST 
for example, will provide far greater sized statistical samples. In- 
deed, we are already beginning to see the benefit of inc reased sam- 
ple si zes for lens statistics from surveys such as SDSS dOguri et"ai] 
l2008h . 

In this work we have demonstrated that usable cosmological 
constraints can be placed on the S7m - f^A plane with 100 time 
delay systems or more (c.f. ~ 20 presently known time delay sys- 
tems), a sample size that is likely to become available within several 
years. When samples grow to several hundreds in the next decade, 
strong lensing has the potential to constrain such cosmological pa- 
rameters to levels approximately equal to that of current CMB ex- 
periments, although of course these will have also increased in pre- 
cision with future measurements. We note that analyses from lens 
statistics produce orthogo nal constraints to those that are derive d 
through time-delays (e.g. [Mitchell etar]|2005l: lOguri et al.ll2008l) . 
revealing the possibility of dual constraints from the same statisti- 



cal sample. This holds further potential in increasing the attainable 
precision for cosmology from strong lensing. 

Currently the best estimate for the Hub ble constant is from th e 
HST Key Project at 72 ± 8 km s~^ Mpc"^ dFreedman et alj200lh . 
This still gives a rather broad 2a range of 56-88 kms"^ Mpc "\ 
a relatively large uncertainty for one of the most fundamental cos- 
mological parameters that governs the length and time-scale of our 
Universe. With the addition of new time delay systems from future 
surveys, lensing has the potential of providing a very high precision 
constraint on Hq. As discussed previously, many of the existing 
lensing degeneracies, e.g. between the Ho and the density profile 
slope, should become less of an issue with large samples and it 
is likely that on average the slope will converge in such a sample 
(i.e ^ 100) to a value which, on current indications, is around 77 = 
2. Indeed, we have seen in this study that larger samples hold the 
promise of reducing the uncertainty of the Hubble constant value 
to levels as low as 1-2% uncertainty with the addition of a fur- 
ther 80 systems to the existing time delay sample. In this study we 
have assumed a mean galaxy density profile that does not evolve 
with redshift, based on the current best estimate from the SLACS 
study, although there still remains an uncertainty in this matter. The 
SLACS study now has over 80 lens systems and dynamical mass es- 
timates of higher redshift systems are being determined - this may 
provide a better empirical understanding of mass slope evolution. 
Large future surveys should be able to resolve this question. 

With upcoming missions such as Planck and the SKA prepar- 
ing to yield improved constraints on many cosmological parame- 
ters it is obvious to question why constraints from lensing should 
be pursued. Lensing as a tool for cosmology is interesting for a 
key reason; one of the most vital questions in cosmology today is 
whether a dark energy component exists. At present only the super- 
novae data set is able to address that question (CMB experiments 
can constrain the total density fim + 0,a, while clusters constrain 
S7m)- A significant result from lens statistics so far is the unique 
and independent evidence it provid es for Q,a > 0, wit hout any 
other cosmological assumptions (e.g. iMitchell et al.l2005l) . The ex- 
act physics that allows Type la supemovae to be used as cosmologi- 
cal probes is still an area of uncertainty, and as such the information 
and confirmation obtained from lensing is highly valuable. Indeed, 
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it is tlie complimentary nature of constraints from strong le nsing 
tliat could prove to be valuable, as discussed in lLindeiH2004h . 

In addition, other interesting possibilities exist to do cosmog- 
raphy with larger strong lensing samples, such as using samples of 
double Einstein rings , as recently discovered in the SLACS data set 
jOavazzi et al. or by co mbining strong lensing and stella r 

dynamics data as suggested bv lGrillo. Lombardi & Bertinl ( |2008|) . 
While here we have focused on w = —1, strong lensing time de- 
lays have bee n proposed to test for quintessence cosmologies (e.g. 
lLindeJl2004h since the equation of state of the quintessence com- 
ponent and its evolution influence the value of the Hubble con- 
stant. Although the differences between evolving and non-evolving 
cosmologies are relatively small, and an order of magnitude more 
lenses than currently available would be needed to ac curately probe 
any quintessence component fcewis & Ibatal l2004h . this will be 
possible in the near future. 

We have seen how lensing time delays have the potential to 
constrain the Qm - J^a plane. However, this analysis depended 
in part upon the assumption of Gaussian distributions for the bulk 
properties of sources and lenses in future lens survey samples. For 
some of the parameters this only served as a necessary approxima- 
tion. Indeed, the nature of the distribution of density slopes, for ex- 
ample, could well have pronounced tails in the "super-isothermal" 
region. While it is true that the average profile of the ensemble may 
be isothermal, a number of systems exhibiting super-isothermality 
could result in non-Gaussian distributions. Similarly, features such 
as ellipticity distributions are not particularly well defined for lens 
galaxies at present. Future work on trying to extract cosmological 
information from large lens time delay samples would do well to 
explore the effects of more realistic, non-Gaussian distributions. 

In this study we have demonstrated what can be expected in 
the way of constraints through the use of a commonly adopted 
elliptical model approach. Different methods such a s the non- 
param etric approach to modelling lens systems (e.g. I Saha et all 
20061), or in deed statistical methods involved in Ho derivations (e.g. 
Oguril200 7) could lead to greater accuracy in the constraints of pa- 
rameters studied here, or similar constraints with smaller sample 
sizes. Whichever method is used in the approach, the new data sam- 
ples from future surveys promise to yield important results from 
strong lensing systems and as such are eagerly anticipated. 
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